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Abstract 

Genetic diversity of 60 Hevea genotypes, consisting of Asiatic, Amazonian, African and lAC clones, and pertaining to 
the genetic breeding program of the Agronomic Institute (lAC), Brazil, was estimated. Analyses were based on 
phenotypic multivariate parameters and microsatellites. Five agronomic descriptors were employed in multivariate 
procedures, such as Standard Euclidian Distance, Tocher clustering and principal component analysis. Genetic vari- 
ability among the genotypes was estimated with 68 selected polymorphic SSRs, by way of Modified Rogers Genetic 
Distance and UPGMA clustering. Structure software in a Bayesian approach was used in discriminating among 
groups. Genetic diversity was estimated through Nei's statistics. The genotypes were clustered into 12 groups ac- 
cording to the Tocher method, while the molecular analysis identified six groups. In the phenotypic and microsatellite 
analyses, the Amazonian and lAC genotypes were distributed in several groups, whereas the Asiatic were in only a 
few. Observed heterozygosity ranged from 0.05 to 0.96. Both high total diversity (H^ = 0.58) and high gene differenti- 
ation (G3,, = 0.61 ) were observed, and indicated high genetic variation among the 60 genotypes, which may be useful 
for breeding programs. The analyzed agronomic parameters and SSRs markers were effective in assessing genetic 
diversity among Hevea genotypes, besides proving to be useful for characterizing genetic variability. 
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Introduction 

The Hevea genus belongs to the Euphorbiaceae fam- 
ily and comprises 1 1 species native to the Amazon region 
(Pires et al, 2002). Hevea hrasiliensis (Willd. ex Adr. de 
Juss.) Muell-Arg. is the only cultivated species and the 
main source of natural rubber. 

Until about 1913, Brazil was the major producer of 
natural rubber, which was obtained from wild rubber trees 
growing in the rain forest of the Amazon basin. However, 
with the introduction of the Wickham material in 1876, 
Southeast Asia has gradually becoine the major producer of 
natural rubber accounting for more than 90% of the total 
production worldwide. There are approxiinately 7 to 8 inil- 
lion hectares of rabber plantations in the rubber areas of 
Asia and Africa. Genetic improvement through mass selec- 
tion and modified recurrent selection has resulted in the 
production and release of elite clones, especially from Ma- 
laysia, over the past 60 years (Onokpise, 2004). Currently, 
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Hevea brasiliensis is cultivated in several tropical coun- 
tries, most of which have active plant-breeding programs 
(Sedgley and Attanayake, 1988). 

In the past, there were limited numbers of H. 
brasiliensis clones suitable for use as parents in breeding 
programs, most of those available having already been se- 
lected according to phenotypic performance. They were 
crossed in many possible combinations, with posterior se- 
lection of the most promising fainilies and progenies. Now- 
adays, an increased number of potential parents are 
available as a result of substantial breeding efforts and the 
exchange of clones among research institutions. Conse- 
quently, a wide range of crosses can now be attempted, this 
requiring additional resources for the effective exploitation 
and wise choice of parental clones. 

Estimates of genetic divergence, through multivariate 
analysis of both agronomic characters and molecular mark- 
ers, should provide valuable data for parent-choice in 
breeding programs. Multivariate analysis based on pheno- 
typic data has been used to assess genetic diversity of rub- 
ber tree (Paiva, 1994; Omokhafe and Alika, 2003), as has 
also occurred with many other plant species, such as the as- 
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sai palm (Oliveira et al., 2007), coffee (Fonseca et al., 
2006) and bean (Chiorato et al., 2007). More recently, mo- 
lecular markers have proved to be useful in estimating ge- 
netic diversity in a wide range of species and populations. 
Among molecular markers, microsatellites or SSRs (Sim- 
ple Sequence Repeats) have received special attention. 
These, besides being codominant and multi-allelic, are 
widely distributed throughout genomes, and thus can be 
highly polymorphic (Chin et al., 1996). Of particular inter- 
est to geneticists and breeders, the SSR markers have been 
successfully used to infer about genetics, pedigree, phylog- 
eny, and/or identity of various traits and/or germplasm ac- 
cessions (McCouch et al., 2001). SSR markers have been 
used to determine genetic diversity in several species, in- 
cluding maize (Laborda et al., 2005), rice (Kwon et al., 
2002), common beans (Benchimol et al., 2007) and rubber 
trees (Lekawipat et al., 2003; Feng et al., 2009). 

The present study reports the suitability of H. 
brasiliensis microsatellite markers, developed from the 
GenBank database, for evaluation of genetic diversity in 
rubber tree clones. Furthermore, the estimates of molecular 
genetic divergence were compared with multivariate phe- 
notypic analysis with the objective of exploring the feasi- 
bility of using SSRs for identifying superior crosses in 
breeding programs. 

Materials and Methods 

Sixty Hevea genotypes (Table 1) from the Rubber 
Tree Program of the Agronomic Institute (lAC, Campinas, 
SP, Brazil) were chosen at advanced evaluation phases 
within genetic breeding programs. The selected genotypes 
consisted of Asiatic, African, Amazonian and lAC clones. 
Several of the Asiatic genotypes were derived from the 
Wickham collection originally introduced into Asia in 
1876, and which are known as Wickham clones. The Ama- 
zonian clones were derived from selection and crossings 
carried out in Brazil by Ford and The North Agronomic In- 
stitute. They are the result of crossings among Amazonian 
and highly productive Asiatic genotypes, with the excep- 
tion of the RO 45 clone, which was derived from a native 
plantation exploited for rubber extraction in the state of 
Rondonia (Brazil). The I AC clones resulted from con- 
trolled crossings and open pollinations performed in this re- 
search institute. 

Phenotypical multivariate statistical analysis 

Average values of five agronomical descriptors, each 
based on three replicates, were subjected to multivariate 
analysis. They comprised average of seven years of girth 
growth increment at juvenile immature phase before tap- 
ping, average of three years of girth growth increment in 
adult trees on tapping, average of three years of dry rubber 
yielding, virgin bark thickness in opened panel tapping; and 
the total number of latex vessel rings. These data were col- 



lected over a period of ten years. Measurements were taken 
as described by Gon9alves et al. (2006). 

Multivariate procedures consisted of Standard Eu- 
clidian Distance (SED), Tocher Clustering and Principal 
Component Analysis (PCA). The contribution of each vari- 
able to genetic divergence was calculated by the criteria of 
Singh (198 1). Statistical analyses were performed using the 
Genes software (Cruz, 2006). 

SSR development and characterization 

Total genomic DNA samples were extracted from 
powdery lyophilized leaf tissues using the 2% CTAB 
method (Hoisington et al. , 1 994) with few modifications. A 
total of 470 reads from GenBank were evaluated in the de- 
velopment and characterization of Hevea microsatellites. 
Redundancies were identified using BLASTN software 
search utilities in GenBank (Ahschul et al, 1990). The SSR 
motifs in the sequences were identified, counted and local- 
ized by using SSRIT (Simple Sequence Repeat Identifica- 
tion Tool_ software. A total of 80 primer pairs (Table SI) 
were developed using Primer Select software from the 
Lasergene program (DNASTar, Inc.). 

PCR amplifications were carried out in a 25 |iL vol- 
ume containing 100 ng of DNA, 1 U Taq DNA polymerase, 
1.5 mM MgCb, 200 |aM of a total dNTP mixture and 
0.8 |iM of each forward and reverse primer. Each SSR was 
characterized on a gradient amplification profile, by vary- 
ing the annealing temperature {Ta) at a difference of up to 
10 °C. After an initial denaturing step of 1 min at 94 °C, the 
PCR amplification was performed in 30 cycles of 1 min at 
94 °C, 1 min at the specific Ta and 1 min at 72 °C, followed 
by a final extension at 72 °C for 5 min and then kept at 
15 °C. Alternatively, some SSRs could only be amplified 
by Touchdown. Amplification products were resolved on 
6% (w/v) denaturing polyacrylamide gels and silver 
stained, according to Creste et al. (2001). 

Polymorphism analysis, genetic distances and SSR 
clustering 

Data on the presence (1) or absence (0) of SSR bands 
were transformed into genotypic data in order to identify 
loci and alleles. The Polymorphism Infonnation Content 
(PIC) value for each locus was calculated using the PIC for- 
mula = "^^f^ fj ^'^^ which n is the number 
/= 1 ,= 1 1 

of alleles; and /J and fj are the frequencies of the * and y* al- 
leles, respectively (Botstein et al. 1980). 

Genetic distances were calculated by using Modified 
Rogers Genetic Distance (MRD) according to Goodman 
and Stuber (1983). A genetic distance matrix was estimated 
using TFPGA software. Cluster analyses were performed 
using UPGMA with the NTSYS-pc computer package ver- 
sion 2.02E. Clustering stability was tested by the Bootstrap 
procedure based on 10,000 re-sampling using the BooD 
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Table 1 - Sixty rabber tree genotypes selected in the breeding program of 





Clone 


Genealogy 


1 


lAC 300 


RRIM 605(Tjir 1 x PB 49) x AVROS 
353(AVROS 164 X AVROS 160) 


2 


lAC 301 


RRIM 501(Pil A 44 X Lun N) x AVROS 


3 


lAC 302 


RRIM 501 (Pil A 44 X Lun N) x AVROS 

J J V IW 1 OH- A /\ V XV w O i DU I 


4 


lAC 303 


RRIM 51 l(Pil A 44 X Pil B 16) X AVROS 
1 1 R r A VRO*\ 91 4 V A VRn*s 7'\f^^ 

I. ^ 1 o y rt. V l\ V./ O ^ i T A rt. V I\ V y O jL^yJ f 


D 


T AM AAQ'i 


Fv AAl 1 i'F A^ll V PR v Tiir 1 
r A ^HJ. I yr ■h-j / j a rri oo ) x ijir i 




IAN 3 1 93 


Fv S 1 f\(V 4.S49^'' y AVROS '\f^'^^ y PR 

1 A J 1 Jt^Z. a rt. V rVv^kj J U J ^ A JT JJ OU 


/ 




AVRn*s 40 V RRTM SflQi'Pil A 44 y T ^^r\^ 


Q 

o 


lAC 307 


A VR n*s 1 9 S f A VR O'i; 9 1 4 V A VR n*s ^ 1 7"! 

xPR 107 


9 


lAC 318 


RRIM 600 (Tjir 1 x PB 49) x Fx 3899 (F 
4542<'>x AVROS 363) 


10 


Fx 3899 


F4542<"x AVROS 363 


11 


IAN 6323 


Tjir 1 X Fx 3810 (F 4542"'x AVROS 363) 


12 


lAC 328 


RRIM 600 (Tjir 1 x PB 86) x PR 107 


13 


lAC 329 


GT711 xTjir 16 


14 


lAC 330 


RRIM 600(Tjir 1 x PB 86) x GT 71 1 


15 


lAC 331 


RRIM 600(Tjir 1 x PB 86) x AVROS 1328 
(AVROS 214 X AVROS 37) 


16 


lAC 333 


C 228 X Tjir 16 


17 


lAC 332 


GT 71 1 X RRIM 600 (Tjir 1 x PB 86) 


18 


lAC 400 


GT 71 1 X RRIM 600 (Tjir 1 x PB 86). 


19 


lAC 401 


RRIM 600 (Tjir 1 x PB 86) ill. 


20 


lAC 402 


GT711 ill 


21 


lAC 403 


GT711 ill. 


22 


lAC 404 


PB 5/63 X AVROS 636 


23 


lAC 405 


Tjir 1 X RRIM 323. 


24 


lAC 406 


IAN 873 (PB 86 x FA 1717) x RRIM 600 
(Tjir 1 X PB 86) 


25 


lAC 407 


RRIM 600 (Tjir 1 x PB 86) ill. 


26 


lAC 408 


RRIM 513 (Pil B 16 X Pil A 44) ill. 


27 


lAC 409 


Fx 2784 (F 4542*" x AVROS 363)111. 


28 


IAC410 


PB 86xPB 235. 


29 


IAC411 


GT711 ill. 


30 


IAC413 


IAN 873 (PB 86 x FA 1717) ill. 


31 


IAC414 


lAC 126 (Fx 25 (F 351 x AVROS 49) x Tjir 
1) ill. 


32 


IAC415 


AVROS 363 ill. 


33 


lAC 412 


IAN 873 (PB 86 x FA 1717) x GT 71 1 


34 


IAC417 


RRIM 600 (Tjir 1 x PB 86) ill. 



program (Coelho, 2002). The cophenetic coefficients be- 
tween the genetic distance matrix and the dendrogram de- 
rived matrix were performed using the NTSYS-pc 
computer package. The significance of cophenetic correla- 
tions was tested by applying Mantel correspondence analy- 
sis. Principal Coordinate Analysis (PCO, Gower 1966) was 
performed using MRD distance matrix. Genetic diversity 
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Instituto Agronomico (lAC) and their respective genealogy. 





Clone 


Genealogy 


35 


lAC 418 


RRIM fton tTWr 1 Y PR RA^ ill 
JVIVllVi ouu ^ iJAi i X t:D oOf ill. 


36 


lAC 420 


TAM 87^ ^PR SA V F A 1717^111 


J 1 


lAC 42 1 


TAP 1 '^7 FFy SnS^F 4^49*'* v AVROS '^f\W v 
Fx 25 (F351 X AVROS 49) ill. 


38 


lAC 422 


RRIM 513 (Pil B 16 X Pil A 44) ill. 


39 


lAC 423 


lAC 90 [RRIM 507(Pil B 84 x Pil A 44) x Fx 
25(F 351 X AVROS 49) ill. 


40 


lAC 424 


RRIM 600 (Tjir 1 x PB 86) ill. 


41 


lAC 425 


RRIM 600 (Tjir 1 x PB 86) ill. 


42 


IAN 3156 


Fx 516 (F 4542'" x AVROS 363) x PB 86 


43 


R0 45 


Primary clone 


44 


lAC 40 


RRIM 608 (AVROS 33 x Tjir 1) x AVROS 
1279 (AVROS 156 x AVROS 374) 


45 


RRIM 701 


44/553 X RRIM 501 (Pil A 44 x Lun N) 


46 


PB235 


PB 5/51 xPB S/78 


47 


GT 1 


Pi'iTnnn/ r*lr»nf* 
JTlllllclly UlUllC 


48 


IAN 873 


PB 86xFA 1717 


49 


RRIM 600 


Tjir 1 X PB 86 


50 


IRCA 130 


PB 5/51 XIR 22 


51 


lAC 15 


RRIM 504 (Pil A 44 X Lun N) x RRIM 600 
(Tjir 1 X PB 86) 


52 


1AC35 


Fx 25 (F 35 1 X AVROS 49) x RRIM 600 
(Tjir 1 X PB 49) 


53 


1AC41 


RRIM 608 (AVROS 33 x Tjirl) x AVROS 
1279 (AVROS 256 x AVROS 374) 


54 


lAC 44 


IAN 2325[PB 86 x Fx 3933 (F4542*" x 

AVROS "^f^^M V AVROS 1 '^7R(' AVROS 7 M 
r\ V JVW^ JO J A rt V ivWl5 i JZO^rt V IS\JJ Z 1 4 

X AVROS 3170] 


55 


PR 255 


Tjir 1 XPR 107 


56 


PR 261 


Tjir 1 XPR 107 


57 


PB217 


PB 5/51 XPB 69 


58 


IAN 3703 


Fx 4371 [Fx 3472(F 4542'" x PB 86) x PB 
86] X PB 86 


59 


IRCA 1 1 1 


PB 5/51 X RRIM 600 (Tjir 1 x PB 49) 


60 


PB 28/59 


PBIG seedling 



' ' Primary clone of Hevea benthamiana. 

Amazonian clones (F = Ford, FA = Ford Acre, Fx = Ford crossbred, 
IAN = Instituto Agronomico do Norte); Clones from the State of Sao 
Paulo (lAC = Instituto Agronomico); Asiatic clones (AVROS = Al- 
gemene Vemeiging Rubber planters Oostkust Sumatra, Indonesia; 
GT = Godang Tapen, Indonesia; PB = Prang Besar, Malaysia; PR = Pro- 
efstation voor rubber, Indonesia; Pil = Pilmoor, Malaysia; RRIM = Rubber 
Research Institute of Malaysia, Malaysia; Tjir = Tjirandji, Indonesia); Af- 
rican clones (IRCA = Institute de Recherches sur le Caoutchouc, Ivory 
Coast). 

among genotypes was estimated by way of Nei statistics us- 
ing FSTAT Software. 

The SED and MRD dissimilarity matrices were corre- 
lated using the Genes software (Cruz, 2006). Both the t- and 
Mantel tests were employed with 10,000 simulations to 
attribute significance values to the data. Intra- and inter- 
group correlations were performed using pair-wise genetic 
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distances within and among groups, separated according to 
the group association pattern observed in the dendrogram 
(Figure 1). 

The Bayesian approach of Pritchard et al. (2000) im- 
plemented by Structure software 2.2, was utilized alterna- 
tively to infer clustering. The number of clusters was 
defined from K=^toK = 20, and ten runs of each K were 
conducted using the admixture model and correlated allele 
frequencies, a 200,000 bum-in period and 500,000 MCMC. 
Ad hoc statistics was related to rate changes in the log prob- 
ability of data according to the number oi Ks proposed by 
Evano et al. (2005), with AA^being used as a predictor of the 
ideal number of clusters. In addition, the ideal number 
according to Pritchard and Wen (2004) was used as the cri- 
terion for defining the number of groups (k). The most 
trustworthy value was estimated based on the lowest nega- 
tive number of Ln (the log-likelihood of the data) and the 
lowest standard deviation found during statistical analysis. 



Group I 



Group! 



Group 3 

Group 4 
Group 5 
Group 6 



1- IAC 300 

2- IAC 301 

3- lAC 302 
■I-IAC 3071 
9.|AC3l8r 

IMAN 63:3 

- 27-IAC-40I) 
■ 1 2-1 AC 328 

13- 1 AC 329 
17-1AC332 

19- IAC401 

is-iAC-mo 

2 1- 1 AC 403 

15- I.\C331 

16- IAC333 

25- IAC407 
.M.IAC4I7 

14- I.AC330 

20- IAC 402 

22- IAC 404 

23- l.« 40.1 
28.IAC410 

26- l.\C 40S 

24- I.AC 406 
33-1.AC41: 
48-lAN 87_! 
32-IAC415 

29- 1AC41 1 

30- IAC413 

35- l.\C418 
37-1 AC 421 

39- IAC 423 
4l-l.^C42.'i 

40- IAC 424 

54- IAC' 44 

36- IAC 420 
.';').1RCA 1 1 1 

46-PI! 235 
hO-Pli :s 5') 

4- t.lAC 40 
41-1A141 

49-RRlM 600 
56-Pr261 
50.1RCA 130 
.Il-IAC I.S 

- 52-l,AC 3.'i 
r 4-IAC 303 

7-IAC 306 
5-IAN 4493 
(i-IAN 3193 
. 10-Fx 3899 
42-lAN 3156 
43-R() 45 
. IAN 3705, 

■45-RRlM 701 
47-GT 1 

55- PR 255 
t 57.PB2I7 

31- IAC4I4 
3S-IAC 422 



Modified Roger's genetic distance 

43.70''. 



] — j55.ir<i 




r 



— r- 

0.36 



1 

0.46 

Coefficient 



— r- 

11.56 



—I 

0.67 



Figure 1 - UPGMA cluster analysis of Modified Rogers Genetic Dis- 
tances based on data from 68 SSRs, used in the evaluation of the 60 rub- 
ber-tree genotypes. Bootstrap node suppoil, represented in percentages, 
shows clustering stability. Nuinbers (%) on the branches correspond to 
bootstrap values above 40% (10,000 replications). (Cophenetic va- 
lue = 0.78). 



Results 

Phenotypic analysis 

The 15 most divergent genotype pairs identified by 
the SED matrix are listed in Table 2. Hevea benthamiana 
was a common ancestor for seven of the 15. The most di- 
vergent genotype pair was lAC 318 - PB 235. 

Through Tocher analysis, the 60 rubber-tree geno- 
types were clustered into 12 groups (Table 3). Among 
these, the eight Amazon genotypes (IAN 6323, IAN 3 156, 
RO 45, IAN 4493, IAN 3193, IAN 873, IAN 3703 and Fx 
3899) were distributed into five groups (I, II, VI, VII, VIII), 
and the 42 lAC genotypes into nine (I, III, IV, V, VI, VII, X, 
XI, XII). These results indicate high genetic diversity in the 
lAC and Amazon genotypes. Diversity in the eight Asiatic 
genotypes (RRIM 701, GT 1, PR 255, PB 217, RRIM 600, 
PR 261, PB 28/59 and PB 235) was low, with the majority, 
except for PB 235, PB 217 and RRIM 701, being allocated 
to group 1 . The two African genotypes (IRCA 111 and 
IRCA 130) were clustered into one and the same group 
(IV). In Group 1, 56% of the genotypes proved to have ei- 
ther of the Indonesian clones, Tjir 1 or GT 71 1, in their an- 
cestry. Notwithstanding, ancestry was not considered to be 
a suitable criterion for characterizing most of the groups. 

PCA forphenotypic data accounted for 80.88% of the 
total variation in the first three principal components. The 
average dry rubber yield was the variable that contributed 



Table 2 - Fifteen pairs of the most divergent genotypes according to Stan- 
dard Euclidian Distance (SED), estimated for 60 rubber-tree genotypes 
and considering five agronomic descriptors. 



Order 


SED 


Genotype pairs 


1° 


6.41 


IAC318-PB 235 


2° 


6.16 


lAC 40 -RRIM 701 


3° 


6.11 


lAC 318 -IRCA 130 


4° 


6.10 


1AC318-1AC 401 


5° 


6.10 


IAC318-1AC 400 


6° 


6.09 


IAN 3156 -PB 217 


7° 


6.08 


LAC 306 -IRCA 130 


8° 


6.06 


IAN 3156 -RRIM 701 


90 


6.05 


lAC 306 - lAC 406 


10° 


5.99 


lAC 306 - lAC 401 


11° 


5.98 


LAC 40 -PB 217 


12° 


5.97 


LAC 414 -PB 235 


13° 


5.85 


lAC 306 - lAC 400 


14° 


5.80 


lAC 414 -IAN 3156 


15° 


5.76 


lAC 331-IAC 401 



Clones from the State of Sao Paulo (lAC = Instituto Agronomico); Asiatic 
clones (PB = Prang Besar, Malaysia; RRIM = Rubber Research Institute 
of Malaysia, Malaysia), Ainazonian clones (IAN = Instituto Agronomico 
do Norte); Clones from the State of Sao Paulo (lAC = Instituto Agro- 
nomico), African clones (IRCA = Institute de Recherches sur le Caou- 
tchouc, Ivory Coast). 
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Table 3 - Clustering of 60 rubber-tree genotypes by the Tocher method 
and based on dissimilarity estimated by Standard Euclidian Distance 
(SED), using five agronomic descriptors. 



Groups 


Genotypes 


I 


lAC 302, Fx 3899, IAN 6323, lAC 328, lAC 329, lAC 330, 




lAC 333, lAC 402, lAC 403, lAC 404, lAC 407, lAC 408, 




lAC 409, lAC 411, lAC 415, lAC 412, lAC 417, lAC 421, 




^ZZ, HjL^, IJ\\^ Kj L i, i/\i>J o/j, KKilVi OUU, 




lAC 15, LAC 35, lAC 41, lAC 44, PR 255, PR 261, 




IAN 3703, PB 28/59 


II 


IAN 3193, RRIM 701, PB 217 


in 


lAC 332, lAC 413, lAC 414, lAC 418 


rv 


lAC 400, lAC 401, IRCA 130, IRCA 1 1 1 


V 


lAC 405, lAC 406, lAC 410, lAC 420 


VI 


lAC 300, lAC 301, lAC 303, IAN 4493, lAC 306, lAC 307 


VII 


RO 45, lAC 40 


VIII 


IAN 3156 


DC 


PB 235 


X 


1AC423 


XI 


1AC318 


XII 


lAC 33 1 



Amazonian clones (Fx = Ford crossbred, IAN = Instituto Agronomico do 
Norte); Clones from the State of Sao Paulo (lAC = Instituto Agronomico 
de Campinas); Asiatic clones (AVROS = Algemene Vemeiging Rubber 
planters Oostkust Sumatra, Indonesia; GT = Godang Tapen, Indonesia; 
PB = Prang Besar, Malaysia; PR = Proefstation voor rubber, Indonesia; 
RRIM = Rubber Research Institute of Malaysia, Malaysia); African clones 
(IRCA = Institute de Recherches sur le Caoutchouc, Ivory Coast). 

the most in the estimation of the genetic divergence among 
the 60 genotypes. The number of latex vessel rings was the 
least important variable, and so could be discarded. 

Molecular analysis 

Of the 80 characterized SSRs (Table SI) 68 were 
polymorphic and informative. In the SSR IAC-Hv34 geno- 
type, two distinct bands were amplified and considered as 
two distinct loci (Hv34a, Hv34b). Thus, 69 polymorphic 
and informative SSR loci were identified and character- 
ized. The electrophoretic profile obtained with SSR lAC- 
-Hv72 can be observed in Figure SI. Polymorphic 
information content varied from 0. 1 1 to 0.87, with an aver- 
age of 0.57. The mean allele number per locus was 5.88, 
ranging from 2 to 13, the extremes being attributed to SSR 
IAC-Hv36 and SSR IAC-Hv20, respectively. 

Six SSRs (IAC-Hv67, IAC-Hv68, IAC-Hv76, lAC- 
Hv44, IAC-Hv69 and IAC-Hv66) were efficient at ampli- 
fying Hevea pauciflora, and showed transferability in 
relation to Hevea brasiliensis, using the same PGR amplifi- 
cation procedure. 

Observed heterozygosity (Ho) based on Nei's esti- 
mates varied from 0.05 to 0.96, with an average of 0.45 (Ta- 
ble 4). Extreme Hq values were encountered in SSR lAC- 
HV09, IAC-HV66 and IAC-HV76. High total diversity 
(Ht' = 0.58) and high gene differentiation (Gst' = 0.61) were 



observed among all the 60 genotypes. SSR PCO accounted 
for 19.66% of the total variation in the first three axes. 
When selecting the most divergent pairwise distance in the 
Rogers modified genetic distance matrix (Table 5), lAC 
414 appeared in eight of the pairwise distances. Den- 
drogram analyses (Figure 1) revealed six distinct groups. 
The eight Asiatic genotypes (RRIM 70 1 , GT 1 , PR 255, PB 
217, RRIM 600, PR 261, PB 28/59 and PB 235) were dis- 
tributed in only two groups (II, V). The African genotypes 
(IRCA 1 1 1 and IRCA 130) were clustered in group II. On 
the other hand, all the lAC genotypes were distributed 
among four groups (I, II, III and VI) while the eight evalu- 
ated Amazonian genotypes (IAN 873, IAN 6323, IAN 
4493, IAN 3193, Fx 3899, IAN 3156, RO 45 and IAN 
3709) were distributed in groups I, II, III and IV. 

Bootstrap analysis expressed high statistical node 
support for genotypes with shorter distances (Figure 2). 
The cophenetic correlation was r = 0.78 (p < 0.002). Groups 
were clearly distinguished, with several clusters being sup- 
ported by high bootstrap values. Bootstrap analysis and 
cophenetic correlations indicated that SSR dendrogram 
clustering accurately depicted estimated genetic distances 
among rubber-tree genotypes. Group 1 contained all the ge- 
notypes derived from the AVROS clones. Group 2 com- 
prised genotypes with either the GT 711 or RRIM 600 
clone in their ascendancy. All the GT 711 derived geno- 
types were clustered in Group 2, and all the RRIM 600, but 
one (lAC 3 1 8), in Group 1 . The lAC 400 genotypes were 
clustered in Group 2, except for two that were positioned in 
Group 6. The Amazonian genotypes IAN 4493, IAN 3193 
and Fx 3899 were clustered in Group 3. Group 4 included 
the other three Amazonian clones IAN 3156, RO 45 and 
IAN 3703. The Wickham clones RRIM 701, GTl, PR 255 
and PB 217 were gathered in Group 5. The lAC 414 and 
lAC 422 clones, placed in Group 6, were the only ones of 
the lAC 400 series outside Group 2. 

A total of six groups were identified by ISK as being 
the ideal number of groups, as previously proposed by 
Evano et al. (2005), and according to criteria indicated by 
Pritchard and Wen (2004). In an investigation of correspon- 
dence between the dendrogram and structure groups 

800000 T 1 
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3 4 5 6 7 8 9 10 II 12 13 14 15 16 17 18 19 20 
K 

Figure 2 - Representation of the number of ideal clusters identified by 
Structure software according to the methodology of Evano et al. (2005). 
The analysis was based on 68 SSRs utilized in the evaluation of the 60 rub- 
ber-tree genotypes. 



Gouvea et al. 313 



Table 4 - Gene diversity analysis of 60 rubber tree genotypes of the breeding program of tlie Instituto Agronomico (lAC) estimated by SSRs. 





rlo 


rlt 


Ost 




rlo 


rlt 


Ost 


IAC-HV09 


0.05 


0.42 


0.94 


IAC-HV45 


0.20 


0.18 


0.45 


lAC-HVlO 


0.30 


0.88 


0.83 


IAC-HV46 


0.55 


0.56 


0.51 


IAC-HV28 


0.27 


0.58 


0.77 


IAC-HV49 


0.36 


0.55 


0.67 


IAC-HV13 


0.52 


0.45 


0.42 


IAC-HV51 


0.73 


0.70 


0.48 


IAC-HV17 


0.40 


0.77 


0.74 


IAC-HV61 


0.68 


0.81 


0.58 


IAC-HV24 


0.33 


0.68 


0.76 


IAC-HV65 


0.12 


0.11 


0.47 


IAC-HV08 


0.23 


0.41 


0.72 


IAC-HV62 


0.39 


0.59 


0.67 


IAC-HV03 


0.14 


0.55 


0.87 


IAC-HV66 


0.05 


0.46 


0.94 


IAC-HV06 


0.47 


0.76 


0.69 


IAC-HV58 


0.57 


0.52 


0.45 


IAC-HV04 


0.55 


0.73 


0.62 


IAC-HV56 


0.06 


0.36 


0.92 


IAC-HV05 


0.46 


0.75 


0.70 


IAC-HV69 


0.90 


0.79 


0.43 


IAC-HV02 


0.53 


0.78 


0.66 


IAC-HV55 


0.43 


0.66 


0.68 


lAC-HVll 


0.43 


0.64 


0.66 


IAC-HV44 


0.65 


0.65 


0.50 


IAC-HV07 


0.77 


0.83 


0.54 


IAC-HV76 


0.96 


0.82 


0.41 


IAC-HV27 


0.47 


0.52 


0.55 


IAC-HV75 


0.41 


0.70 


0.71 


IAC-HV20 


0.65 


0.84 


0.61 


IAC-HV78 


0.60 


0.50 


0.40 


IAC-HV15 


0.62 


0.80 


0.62 


IAC-HV63 


0.32 


0.57 


0.72 


IAC-HV12 


0.27 


0.29 


0.53 


IAC-HV68 


0.17 


0.34 


0.76 


lAC-HVOl 


0.42 


0.61 


0.66 


IAC-HV67 


0.32 


0.49 


0.68 


IAC-HV16 


0.70 


0.72 


0.51 


IAC-HV79 


0.33 


0.34 


0.52 


IAC-HV23 


0.77 


0.60 


0.36 


IAC-HV70 


0.21 


0.57 


0.82 


IAC-HV25 


0.55 


0.65 


0.58 


IAC-HV50 


0.33 


0.34 


0.52 


IAC-HV29 


0.25 


0.25 


0.48 


IAC-HV57 


0.42 


0.40 


0.48 


IAC-HV22 


0.82 


0.81 


0.49 


IAC-HV80 


0.64 


0.62 


0.48 


IAC-HV30 


0.87 


0.84 


0.49 


IAC-HV53 


0.88 


0.59 


0.26 


IAC-HV18 


0.72 


0.72 


0.51 


IAC-HV74 


0.82 


0.85 


0.52 


IAC-HV14 


0.68 


0.75 


0.55 


IAC-HV75 


0.58 


0.60 


0.52 


IAC-HV31 


0.48 


0.73 


0.67 


IAC-HV52 


0.29 


0.39 


0.63 


IAC-HV32 


0.10 


0.13 


0.60 


IAC-HV73 


0.60 


0.84 


0.64 


IAC-HV40 


0.36 


0.51 


0.65 










IAC-HV33 


0.67 


0.79 


0.58 










IAC-HV35 


0.28 


0.67 


0.79 










IAC-HV42 


0.38 


0.53 


0.64 










IAC-HV47 


0.35 


0.49 


0.65 










IAC-HV34a 


0.90 


0.59 


0.24 










IAC-HV34b 


0.15 


0.14 


0.47 










IAC-HV49 


0.48 


0.65 


0.63 










IAC-HV38 


0.07 


0.16 


0.79 










IAC-HV60 


0.17 


0.56 


0.85 










Overall 


0.45 


0.58 


0.61 











Ho: observed heterozygosity. 

Ht' : total heterozygosity. 

Gst' : co-efficient of gene differentiation. 



(Figure 3), group 1 of the dendrogram corresponded entirely 
to Structure group 1, and included genotypes derived from 
crosses of enhanced clones from the Rubber Research Insti- 
tute of Malaysia (RRIM) and Algemene Vemeiging Rubber 



planters Oostkust Sumatra of Indonesia (AVROS). Group 2 m 
the dendrogram corresponded to Structure groups 2, 3 and 5. 
Dendrogram group 3 corresponded to Structure group 4, and 
was characterized by Amazonian clones. Dendrogram group 
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Figure 3 - Distribution of the 60 rubber-tree genotypes in groups according to Structure analysis (k = 6), based in 68 SSRs used in the evaluation of the 60 
rubber-tree genotypes. The individuals were represented by vertical bars, each color being associated to a different group. Genotype identification is re- 
ferred to in Table 1 . 



Table 5 - Fifteen pairs of the most divergent rubber-tree genotypes ac- 
cording to Modified Rogers Genetic Distance (MRD) estimated among 60 
rubber-tree genotypes evaluated with 68 SSRs. 



Order 


MRD 


Genotype Pairs 




0.74 


1AC 328 -PB217 


2nd 


0.73 


1AC414-IAC41 


3"" 


0.73 


1AC414-PB217 


4"^ 


0.72 


1AC418-RRIM 701 


5* 


0.72 


1AC 330-RRIM 701 




0.72 


Fx 3899 -PB 217 




0.72 


lAC 328 - lAC 414 


8* 


0.71 


Fx 3899 - PR 255 


9* 


0.71 


lAC 408 - lAC 414 


10* 


0.71 


lAC 422 - IRCA 130 


11'" 


0.71 


lAC 401 -IAC414 


12* 


0.70 


1AC 333 -RR1M 701 


13* 


0.70 


1AC331 -IAC414 


14* 


0.70 


IAN 4493 -lAC 414 


15* 


0.70 


1AC414-PR261 



Clones from the State of Sao Paulo (lAC = Instituto Agronomico de 
Campinas); Asiatic clones (PB = Prang Besar, Malaysia; PR = Proefsta- 
tion voor rubber, Indonesia; RRIM = Rubber Research Institute of Malay- 
sia, Malaysia); Amazonian clones (Fx = Ford crossbred; IAN = Instituto 
Agronomico do Norte); African clones (IRCA = Institute de Recherches 
sur le Caoutchouc, Ivory Coast). 

5 and 6 corresponded to Structure group 6, and included four 
clones of the Wickham collection. 

An interesting clustering aspect, as revealed by the 
structure program, was the distribution of the Amazonian 
and lAC genotypes into several groups, viz., I, III, IV, V 
and I, II, III, IV, V, VI, respectively, whereas the Asiatic 
genotypes (RRIM 701, GT 1, PR 255, PB 217, RRIM 600, 
PR 261, PB 28/59, PB 235) were distributed into only two 
groups (III, VI). This distribution pattern is in agreement 
with the data obtained through Tocher analysis of pheno- 
typic data and the UPGMA dendrogram based on SSRs. 



Matrix correlation between both kinds of genetic dis- 
tances was significant by t-test and Mantel test (r = 0.13, 
p < 0.01). Genetic distances estimated from phenotypic and 
molecular traits were correlated. Pair- wise distances within 
and among groups were separated in the dendrogram ac- 
cording to the respective group association pattern (Figu- 
re 1). Significant values were found for intra-Group 2 
correlations (r = 0.165, p < 0.01) and for inter-Group 1x3 
correlation (r = 0.565, p < 0.01); inter-Group 1x5 (r = 0.547, 
p < 0.01) and inter-Group 1x6 (r = 0.620, p < 0.05). 

Discussion 

Phenotypic analysis 

Standard Euclidian Distances (SED) detected higher 
divergence between the clone lAC 318 and the Asiatic 
clone PB 235, the latter having been derived from a cross- 
ing between two Malayan clones. This dissimilarity possi- 
bly occurred through lAC 318 has the clone Fx 3899 as 
male parental, which is an interspecific hybrid of H. 
benthamiana x H. brasiliensis (Table 1). Among the 15 
pairs of most divergent genotypes (Table 2), H. 
benthamiana appears as the ancestor in seven pairs. Hy- 
bridization may have several evolutionary consequences, 
these possibly including increased intra-specific genetic di- 
versity (Rieseberg, 1997). 

Total variance (80.8%), as explained by the three 
principal components of the phenotypic data, was expres- 
sive. However, it was less than that described by Paiva 
(1994), consisting of 94.76% and 97.49% data variance in 
the first three components with eight and seven descriptors, 
respectively. The number and nature of variables certainly 
have to be taken into consideration on comparing relative 
final variance. 

IVIolecular analysis 

Polymorphic information content was high for SSR 
loci, and indicated a substantial genetic information content 
in the clones analyzed with microsatellites. Microsatellite 
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marker analysis is very efficient when examining genetic 
diversity (Laborda et al, 2005; Saha et al, 2005). PIC val- 
ues for SSR loci were superior to those observed by Feng et 
al. (2009) when using EST-SSRs to analyze cultivated 
clones in rubber trees. Accordingly, PIC values ranged 
from 0 to 0.684 and averaged 0.383. As expected, EST- 
SSRs have been reported as being less polymorphic than 
genomic SSRs in crop plants due to DNA sequence conser- 
vation in transcribed regions (Scott et al, 2000; Eujayl et 
al, 2001). The mean allele number found in this study was 
5.88, varying from 2 to 13. This result was similar to the 
previously reported 5.92, which varied from 3 to 10, in a set 
of cultivated genotypes (Lekawipat et al. , 2003), and higher 
than the average of 2.47 alleles observed by Feng et al 
(2009), when using EST-SSRs. 

Gene differentiation (Gst' ^ 0.61) was high, indicat- 
ing 61% of total variation to be exploited among the avail- 
able genotypes. A total of 52 SSRs presented high 
coefficients of genetic differentiation (Gst' > 0.50). This 
value was higher than that observed for other open polli- 
nated species belonging to the same family as rubber trees, 
such as cassava (Fregene et al, 2003; Lokko et al, 2006). 
Observed heterozygosity was more variable than that de- 
scribed by Saha et al (2005) in a cultivated rubber tree ge- 
notype when using four SSR markers, and therefore 
considered highly informative. This probably accounts for 
the lower variation observed. 

Genomic transferability of the SSR loci between H. 
brasiliensis and H. pauciflora indicated that these SSRs 
could be useful for studies of synteny within the Hevea ge- 
nus. Saha et al (2005) also observed that SSRs specifically 
developed for H. brasiliensis efficiently amplified H. 
benthamiana and H. spruceana, thereby implying the high 
conservation of flanking microsatellite genomic regions. 
More recently, Feng et al (2009) developed EST-SSRs for 
H. brasiliensis, and observed interspecies transferability by 
amplifying//, spruceana, H. nitida,H. benthamiana and//. 
pauciflora, and intergenus transferability in castor oil 
plants (Ricinus communis L.) and cassava (Manihot 
utilissima). 

What makes the Bayesian approach interesting in the 
study of population genetic structures is the facility in de- 
tection without the need for prior information on individual 
origin (Pritchard et al. , 2000). When considering all the in- 
dividuals and clustering in Structure atK ^ 6, the arrange- 
ment was such as to correspond to dendrogram clustering. 
The coherence in genotype clusters indicated a non-random 
distribution of alleles and their frequencies. In fact, clones 
in most groups were gathered according to ascendancy as 
previously described. The lack of consistency in various 
dendrogram clusters could be associated to low bootstrap 
node support (with 10,000 re-samplings) in some of the ma- 
jor groups. Low bootstrap nodes could be associated with 
the lack of genetic structure, and the incapacity of a clearing 
clustering tendency supported by high cophenetic values. 



Phenotypic and Molecular analysis 

On comparing phenotypic and molecular clustering 
patterns (Table 3; Figure 1), molecular markers have 
proved to be very efficient in group characterization by ge- 
nealogy. In rubber tree studies using molecular markers 
(Varghese et al, 1997; Feng et al, 2009), genealogy has 
also been used as an aid in group characterization, although 
not always with satisfactory results. In this work, illegiti- 
mate genotypes derived from open pollination with the 
same female parental, whereas those legitimate were the re- 
sult of controlled pollination, with both common parental 
belonging to separate groups. According to Varghese et al 
(1997) Hevea being a predominantly cross-pollinated tree 
species, has Fl hybrids fixed vegetatively while clones are 
highly heterozygous. As a result of segregation and inde- 
pendent assortment in these clones, the proportion of mar- 
ker alleles in the Fl hybrid from each parent can vary 
considerably. Thus, in highly heterozygous species with a 
common ancestry, pedigree information may not always re- 
veal the exact nature of genetic relationships. 

In phenotypic and microsatellite analyses, both the 
Amazonian and lAC genotypes were clustered into several 
groups, thereby indicating high genetic diversity among 
these genotypes. On the other hand, Asiatic genotypes were 
distributed in only a few groups, thus indicating low diver- 
sity, with, most certainly, a narrow genetic base (Besse et 
al, 1994; Varghese et al, 1997). Molecular analysis indi- 
cated the close relationship between African and Asiatic 
genotypes, thereby confirming their Asiatic genetic base 
(Besse et al, 1994). The data confirmed a narrow genetic 
base for the Asiatic and African genotypes and high genetic 
variability for the Amazonian genotypes. This higher ge- 
netic variability in wild Amazonian genotypes was ex- 
pected, although it is not always associated to desirable 
agronomic breeding characters. Clement-Demange et al. 
(2001) reported that Amazonian wild genotypes do not al- 
ways contribute desirable traits to rubber tree genetic 
breeding. Even though, Amazonian wild genotypes were 
introduced into rubber tree breeding by crossing, so as to 
widen the genetic basis of Asiatic clones. The Amazonian 
genotypes evaluated in this study arose from prior breeding 
and selection procedures, thereby possibly constituting an 
interesting genetic background to be exploited in rubber 
tree breeding programs. Indeed, these clones showed the 
highest genetic diversity when compared to the other im- 
proved clones analyzed, and appear to be attractive for rub- 
ber breeding. 

The most divergent genotypes identified in SED anal- 
ysis (Table 2) differed from those identified through MRD 
analysis (Table 5). H. benthamiana was an ancestor in 
seven of the 15 most dissimilar SED genotype pairs. More- 
over, the rubber tree clone lAC 414 was the most dissimilar 
MRD genotype, being involved in eight of the 1 5 largest 
distances registered. It should be noticed that lAC 414 has 
the Amazonian Fx 25 clone as an ancestor. Although Asi- 
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atic genotypes ancestry was predominant in the majority of 
the evaluated rubber tree clones, higher genetic divergence 
was observed in those genotypes derived from inter-spe- 
cific crosses performed in the past, as well as those having 
the Fx 25 clone as ancestor. The data indicated the strong 
contribution of these two genotypes to genetic divergence 
as described in this study. Further analysis should be under- 
taken to confirm these findings. 

Although both phenotypic and molecular analysis 
revealed differences in genotype clustering, they shared 
several common aspects, such as high diversity between 
Amazonian and I AC genotypes and low among Asiatic. 
Most likely, the difference in genotype clustering was due 
to agronomic descriptors being associated to gene expres- 
sion and may suffer environmental influence modulating 
the phenotype. In contrast, molecular markers, including 
microsatellites, are mostly neutral and consist of either ex- 
pressed DNA sequences, or non-expressed genomic re- 
gions such as introns or regulatory sites. In agreement with 
Moser and Lee (1994), as a result of the complex nature of 
the relationship between phenotypic and genotypic varia- 
tion, genotypes that are phenotypically different may vary 
at only a few loci, and those very similar in appearance or 
performance may be quite dissimilar genetically. In addi- 
tion and as reported by Grivet and Noyer (2003), the rela- 
tionship between neutral polymorphism revealed by mark- 
ers and the polymorphism of useful morpho-agronomic 
characters is not clear. Markers revealed sequence similar- 
ity between individuals in a sample of a locus. Morpho- 
agronomic characters measure resemblances between indi- 
viduals based on variables whose level of expression de- 
pends on the number of potentially epistatic genes. 

Despite being significant, matrix correlation for all 
genetic distances of different marker systems exhibited 
poor association, which was also observed by Roldan-Ruiz 
et al. (2001). Taking into account the molecular dendro- 
gram clustering and correlated genotypes, significant asso- 
ciations were detected between intergroup genetic 
distances (phenotypic and SSR-based GDs) in Groups 1x3; 
1x5 and 1x6. Currently, associations between molecular 
and phenotypic data tend to be stronger in crosses between 
genotypes of similar pedigrees (Smith et a!., 1990). Cer- 
tainly, correlations between phenotypic-based and molecu- 
lar-based distances will be improved as additional probes 
or marker loci are employed in analysis (Moser and Lee, 
1994). 

The genotypes analyzed in the present work had been 
previously selected and evaluated in breeding programs, 
and had demonstrated good performance. Thus, prior 
screening of the most divergent genotype pairs identified 
through both methods is suggested for evaluation of the rel- 
ative agronomic performance of their hybrids. In a tradi- 
tional breeding program, thousands of crosses are normally 
performed and evaluated in experimental designs. Accord- 
ing to the data described in this work, SSR-based genetic 



distances could be useful in selecting superior crosses be- 
tween rubber tree clones derived from a population with a 
broad genetic base. Hence, the application of SSR markers 
in rubber-tree breeding could be instrumental in reducing 
the number of single-cross hybrids to be evaluated. SSRs 
are easily assayed by the Polymerase Chain Reaction, and 
have proved to constitute a potent tool for characterizing 
genetic diversity. Furthennore, in perennial plant species, 
such as rubber trees, this molecular marker technology has 
the additional advantage of shortening breeding time by al- 
lowing for the screening of seedlings and juvenile plants. 
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